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ABSTRACT 

Stellar mass black hole X-ray binaries exhibit X-ray spectral states which also have 
distinct and characteristic temporal properties. These states are believed to correspond 
to different accretion disk geometries. We present analysis of two long XMM-Newton 
observations of the Ultra-Luminous X-ray source (ULX) NGC 1313 X-l, which reveal 
that the system was in two different spectral states. While spectral variations have 
been observed in this source before, this data provides clear evidence that the spectral 
states also have distinct temporal properties. With a count rate of ~ 1.5 counts/s and 
a fractional variability amplitude of ~ 15%, the ULX was in a high flux and strongly 
variable state in March 2006. In October 2006, the count rate of the ULX had reduced 
by a factor of ~ 2 and the spectral shape was distinctly different with the presence of 
a soft component. No strong variability was detected during this low flux state with an 
upper limit on the amplitude < 3%. Moreover, the spectral properties of the two states 
implies that the accretion disk geometry was different for them. The low flux state 
is consistent with a model where a standard accretion disk is truncated at a radius 
of ~ 17 Schwarzschild radius around a ~ 200M Q black hole. The inner hot region 
Comptonizes photons from the outer disk to give the primary spectral component. 
The spectrum of the high flux state is not compatible with such a geometry. Instead, 
it is consistent with a model where a hot corona covers a cold accretion disk and 
Comptonizes the disk photons. The variability as a function of energy is also shown 
to be consistent with the corona model. Despite these broad analogies with Galactic 
black hole systems, the spectral nature of the ULX is distinct in having a colder 
Comptonizing temperature (~ 2 keV) and higher optical depth (~ 15) than what is 
observed for the Galactic ones. 

Key words: Accretion, Accretion Disks, Black Hole Physics, X-Rays: Binaries, X- 
Rays: Galaxies, X-rays: individual: NGC 1313 X-l 



1 INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are off-nuclear, 
compact X-ray sources with luminosities exceeding ~ 
10 39 erg s _1 . The nature of ULXs continues to be an enigma, 
since their isotropic high energy output surpasses the Ed- 
dington limit of even the most massive stellar mass black 
holes, sometimes by large factors. One individual ULX may 
outshine the rest of the galaxy in the high energy X-ray 
band. Several models have been proposed to explain the 
high luminosities of ULXs. They may be powered by "in- 
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termediate mass blac k hole (IMBH)" with m asses Mbh — 
10 2 - 10 4 Mq (e.g., IColbert fc Milled 120051 ) which bridge 
the gap between stellar mass black holes in X-ray binaries 
(XRBs) and super-massive black holes in active galactic nu- 
clei. On t he other hand, th ey may be XRBs with anisotropic 
emission dKing et al.ll200lh or beamed XRBs with relativis- 
tic jets directly pointing towards us, i.e., sc aled down ver- 
sions of blazars (jMirabel fc Rodrigued ll999h Alternatively, 
they may be st ellar mass black holes with s uper-Eddington 
accretion rates (jBegelma n 2002; King 2008). 

In the absence of a direct measurement of the black hole 
mass in ULX, indirect evidence is provided by the spectral 
and temporal properties of these sources. Qualitative and 
quantitative similarities (and differences) between them and 
Galactic X-ray binaries harboring stellar mass black holes 
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can provide insight into their nature. The disco very of orbital 
modulations from the brigh t ULX M 82 X-l |Kaaret et al.l 
l200d : iKaaret fc Fensl l2007h and possibly from a ULX in 
IC 342 (|Sugiho et al.ll200ll ) implies that they too are binary 
systems. 

Galactic X-ray binaries exhibit spectral states (e.g. see 
Zd ziarski fc Gierlihskill2004l . for a review). Detailed spectral 
modeling of these systems reveal that the accretion disks 
have different geometries for different states. In the high 
soft state, the spectra are dominated by emission from a 
cold accretion disk which extends probably to the inner most 
stable orbit. A corona or active region above the cold disk 
comptonizes disk photons to produce hard X-rays. In the 
low hard state, the disk is truncated at a larger radius and 
the X-rays originate from an hot inner disk. In the very high 
state, the geometry is probably similar to that of the high 
state, but here the flux from the corona is nearly equal to 
that of the disk. X-ray binaries are variable on a wide range 
of time-scales and exhibit broad band noise as well as narrow 
features called Quasi-periodic oscillations (QPO). Spectral 
state transitions are always associated with distinct changes 
in the rapid temporal behavior of the systems. 

Both Chandra and XMM-Newton observations of ULXs 
have sometimes shown soft X-ray excess emission which 
has been interpreted as the optically thick emission from 
thin accretion disks with temperatures in the range of ~ 
100 — 300 eV, which indicates that they are I MBH accret- 
ing at sub-Eddington (~ 0.1L E dd) rates ( e.g., iMiller et al.l 
l2003l ; lKaaret et al.ll2003l ; iMiller et al.ll2004r i. While this sug- 
gests that ULX may be scaled up versions of X-ray binaries, 
it should be noted that there are important spectral differ- 
ences. Unlike X-ray binaries, a high energy spectral curva- 
ture has been detected in the the high quality X-ray spec- 
tra of a numb e r of b right ULXs dAgrawal fc Misra 2006; 
Stobbart et al.l 120061 ; iGoad et all 120061 ; iDewangan et"ail 



20061 ; lRobertsH2007l ). This spectral curvature is consistent 



with strong, nearly saturated Comptonization. Such strong 
Comptonization models have been invoked to explain the 
very high state of black hole binaries such as XTE 1550-564 
l|Done fc Kubotal 120061 ). However, the temperature derived 
for the electron plasma (~ 2 keV) is significantly smaller 
than that of X-ray binaries (> 20 keV). 

Another striking similarity between ULX and X-ray bi- 
naries is the presence of significant variability power as well 
as featu res in the power density sp e ctra such as QPOs and 
break s (IStrohmaver fc Mushotzkvl 120031; IDewangan et alj 
l2006l l?l; lMucciarelli et al.ll2006l ; IStrohmaver et al.ll2007r i. The 
QPO and break frequencies are found to be smaller than 
those of X-ray binaries which is taken as evidence that ULX 
harbor IMBH. For some sources, like the bright ULX Holm- 
berg II X-l, the v ariability is absent or very weak (< 2%) 
jGoad et alj|2006h . 

Studies of a collection bright ULX in nearby galaxies re- 
veal that for some ULX a power-law spectral fit (similar to 
the low/hard state of X-ray binaries) is adequate while oth- 
ers require a thermal component (similar to the high/soft 
state). The temperature of the thermal component is ei- 
ther ~ 0.2 or ~ 2 keV suggestive o f a further subdivision 
jWinter et al.ll2006l : [Devi et a.l.ll2007h . These studies do not 
indicate whether these spectral differences represent differ- 
ent spectral states or if ULX themselves are a heterogeneous 
class comprising different kinds of sources. Long term spec- 



tral variability, which is interpreted as state transitions sim- 
ilar to X-ray binaries, has been reported for several ULX. 
Repeated Chandra observations of the Antennae have shown 
hardness ratio changes in several ULXs. Transition from a 
high soft state to a l ow hard state has b een observed from 
two ULXs in IC 342 ilKubota et al.ll20o"ll ) and Holmberg IX 
X-l (|La Parola et al.l 200lf ). High/hard to low/soft spectral 
variability has also be en observed from Holmberg II X-l 
(|Dewangan et al.ir2004h . 

However, these spectral variations have not been shown 
to be associated with temporal behavior changes and hence 
their identification as spectral state transitions similar to 
those observed in X-ray binaries remains uncertain. More- 
over, the spectral variations recorded do not necessarily 
imply a change in the accretion disk geometry. Well con- 
strained spectral parameters are needed for such an infer- 
ence. Hence, good high quality data from two long obser- 
vations of a source undergoing such changes are needed to 
ascertain whether ULX undergo true state transitions. 

NGC 1313 is a s piral galaxy and is lo cated at a distance 
of 4.13 ± 0.11 Mpc (IMendez et al.ll2002l l. The galaxy hosts 
three bright and well separated X-ray sources: two ULXs X- 
1 and X-2, and a supernova remnant SN 1978K. NGC 1313 
X-l is one of the first ULX observed by XMM-Newton to 
show soft X-ray exces s emission that le d to the ~ IOOOM0 
IMBH interpretation (|Miller et alj|2003l ). The ULX belongs 
to the low-temperature, high-luminosity class in the lumi- 
nosity versus temper ature (L — kT) diagra m of bright ULXs 
in nearby g alaxies (|Feng fc Kaaretl 120051 ). iFeng fc Kaaret] 
l|2006l . 120071 ) have utilized the first 14 XMM-Newton obser- 
vations and reported long term spectral variations of the 
two ULXs NGC 1313 X-l and X-2. For X-l, they fitted the 
spectra with an absorbed power-law model and reported a 
correlation between the photon index, F and luminosity. An 
additional soft black body component was required for six of 
the observations. The short exposures of these observations 
(in the range of 9 — 41 ks), did not allow for detailed spectral 
and temporal studies. 

Here we present detailed temporal and spectral proper- 
ties of the ULX using two XMM-Newton observations of 123 
and 20 ks . We describe the XMM-Newton observations and 
the spectral and temporal analysis in §2. In §3, the accretion 
disk geometries inferred from the analysis are described. We 
summarize and discuss the results in 54. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed NGC 1313 with XMM-Newton for a 123 ks long 
exposure on 15 October 2006. NGC 1313 was also observed 
on 6 March 2006 with an exposure of 21.8 ks. We used SAS 
8.0.1 with updated calibration to process and filter both 
EPIC-pn and MOS data. The long observation of 15 Octo- 
ber 2006 was partly affected with low amplitude flaring par- 
ticle background towards the end of the exposure. Cleaning 
of the flaring particle background resulted in 89.3 ks data. 
The observation of 6 March 2006 was not affected by par- 
ticle background. For temporal analysis we selected events 
with patterns — 12 and used either continuous exposure 
free of particle background or corrected for the background 
contribution. For spectral analysis, we selected the events 
with pattern — 4 (single and double) for the EPIC-pn and 
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Figure 1. The background corrected EPIC-pn lightcurves of NGC 1313 X-l binned with 200 s in the (a) 0.2 - 10 keV, (b) 0.2 - 2 keV, 
(c) 2 — 10 keV, and (d) the hardness ratio derived from the observations of 6 March 2006 (left panel) and 15 October 2006 (right panel). 
The relative contribution of background emission in the source extraction region is also shown for the 0.2 — 10 keV band. 



— 12 for the MOS and excluded events adjacent to the bad 
pixels with FLAG=0. For the observation of March 2006, 
the source extension overlaps with the EPIC-pn chip gap 
and a large fraction of events have uncertain pattern. These 
events with poor spectral calibration were excluded from the 
spectral analysis but were included in the temporal analysis. 



2.1 Temporal Analysis 

We extracted two sets of lightcurves from each of the obser- 
vations with bin sizes of 200 s and 0.5 s. A circular region 
of 40" centered at the source position was used to extract 
the EPIC-pn source lightcurves. We also extracted back- 
ground lightcurves using four circular regions of radii 40" 
in the source-free areas and close to the source and used 
them to correct the source lightcurves for background con- 
tribution. In Figure [1] we show the background corrected 
lightcurves of NGC 1313 X-l in the 0.2-10 keV, 0.2-2 keV 
and 2 — 10 keV bands. For comparison, we have also plotted 
the background contribution in the source extraction region. 
Evidently background contribution is nearly constant and 
negligible throughout the observations. 

As can be seen in Fig. [Q (left panel), NGC 1313 X-l 
showed strong X-ray variability throughout the observation 
of 6 March 2006. The most remarkable variability in the 
0.2 — 10 keV band occurred at an elapsed time of ~ 13000 s 
when the net count rate decreased by a factor of about two 
in 3000 s. The hardness ratio, calculated as the ratio of count 
rates in the 0.2 — 2 keV and 2 — 10 keV band is also variable 
and correlated with the lightcurve. Thus, at high flux, the 
ULX has a harder spectrum. 

In contrast to the 20 ks observation on 6 March 2006, 
the long observation of 15 October revealed that (i) the ULX 
was in a low flux state, X-ray emission dropped by a factor 



of about ~ 2 , and (ii) the emission was remarkably steady 
throughout the ~ 120 ks exposure. The soft and hard band 
emission and the hardness ratio also remained steady, im- 
plying that the spectrum of the ULX did not change sub- 
stantially (see Fig. [1] right panel). Henceforth, throughout 
this work we refer to the 6 March observation as the high 
state and the 15 October one as the low state. 

To quantify the variability, we have calculated the in- 
trinsic source variability expressed in terms of the frac- 
tional variability amplitude (F var = (<j\ s /x 2 ) 1 ^ 2 , where 
o\s i s the variance after subtracting the contribution ex- 
pected fr om measurement err ors, and x is the mean count 
rate ( see IVaughan et alj 1200 j . and references therein). We 
used the background corrected lightcurves binned with 200 s 
to calculate F var , listed in Table [T] While the measured 
F var is consistent with no intrinsic variability during the 
low state, significant intrinsic variability is seen during the 
high state. We have also derived the power density spec- 
tra using the EPIC-pn lightcurves binned with 5 s. For the 
low state, we used the first ~ 79 ks continuous exposure in 
order to avoid the data gap and the flaring particle back- 
ground seen towards the end of the observation (see Fig. [TJ. 
The light curves were divided into five equal segments and 
the PSDs of each segment was computed. The five individ- 
ual PSDs were averaged and the errors were assigned to 
be the standard deviation at each frequency bin. The re- 
sultant PSDs (after binning) are shown in Figure [2] The 
PSD shows red-noise behavior at low frequencies for the 
hard state, while it is flat during the low one. A constant 
model fitted to the PSD for the low state resulted in a mini- 
mum x 2 — 59.3 for 52 degrees of freedom (dof). The best-fit 
constant power density is 2. 12 ±0.04 which is slightly higher 
than the power density of 2 expected from pure Poisson noise 
arising from photon counting statistics. The marginal higher 



4 G. C. Dewangan et al. 

Table 1. Variability properties of NGC 1313 X-l 
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0.2 - 10 keV 


1.47 ±0.09 


0.016 ±0.004 


14.6 ±0.6 


0.77 ±0.07 


0.026 ± 0.006 


9 n+0.9 
z -°-1.9 


0.2 - 2 keV 


1.06 ±0.08 


0.009 ± 0.003 


11.0 ±0.8 


0.55 ± 0.06 


0.010 ±0.004 
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2-10 keV 


0.41 ± 0.05 


0.006 ± 0.003 


23.7 ± 1.2 


0.22 ± 0.04 


0.013 ±0.004 


< 4.3 



Calculated from background subtracted lightcurve with 200 s bins. 
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Figure 2. Power density spectra of NGC 1313 X-l derived from the EPIC-pn lightcurves obtained on 6 March 2006 and 15 October 
2006 



power density could arise from slight variation in the back- 
ground level (see Fig. [T] right panel). Thus, we measure no 
significant intrinsic variability, consistent with the measured 
Fvar (see Table [TJ . As expected, a constant power density 
model fit to the PSD of the high state resulted in a statis- 
tically unacceptable fit (x 2 /dof — 211.5/38). Addition of a 
power-law component (oc v~ a ) improved the fit significantly 
ix 2 /dof = 47.7/36) with a = -2.25 ± 0.30. Interestingly a 
broken power-law model further improved the fit marginally 
(x 2 /dof = 40.9/34). However, the parameters of the model 
are not well constrained and hence we conclude that the 
data is suggestive of a break, but no concrete inference can 
be made. 

2.2 Spectral Analysis 

We extracted EPIC-pn and MOS spectra from the two ob- 
servations using circular regions with radii of 40" centered 
at the position of NGC 1313 X-l. We also extracted back- 
ground spectra using appropriate nearby circular regions 
free of sources. We created spectral response files using the 
SAS tasks rmfgen and arfgen. Spectral bins were chosen 
such that there was a minimum of 20 counts per spec- 
tral channel for the pn and the MOS data of March 2006 
(high state). The high signal-to-noise EPIC-pn, MOS1 and 
MOS2 data extracted from the deep observation of Octo- 



ber 2006 (low state) were grouped to minimum counts of 
200, 50 and 50, respectively. These spectra were analyzed 
with ISIS 1.4.9. All spectra were fitted in the energy range 
0.3 — 10 keV and the errors on the best-fit spectral param- 
eters are quoted at a 90% confidence level. 

We begin with the spectral analysis of the low state. 
A simple absorbed power law (PL) model was fitted indi- 
vidually to the pn and MOS spectra. Examination of the 
individual fits showed the photon indices and absorption 
columns obtained for three data sets were similar within 
errors and the fits had similar residuals. Therefore, we 
fitted the pn and MOS data jointly and used an over- 
all normalization constant to account for possible differ- 
ences in source extraction areas or calibration uncertain- 
ties. The simple power-law model resulted in an statisti- 
cally unacceptable fit (x 2 /dof = 1484.0/674). The ratio 
of the pn±MOS data and the power-law model is shown 
in Fig. ^bottom panel). This plot clearly shows a broad 
hump or a cutoff near 6 keV and a soft X-ray hump or 
excess emission in the 0.6 — 1.2 keV band. Similar spectra 
have been reported earlier for several ULXs (jFeng fc Kaareti 
120051 ; IStobbart et~ai1l2006l ; iDewangan et al.ll2006l ). Addition 
of a multicolor disk blackbody (MCD) component to the PL 
model improved the fit significantly (x 2 /dof = 807/672). 
Replacing the PL with a cutoff PL further improved the 
fit (x 2 /dof = 717.3/671). This kind of a high energy cut- 



State transition of the ULX NGC1313 X-l 5 




Energy (keV) 

Figure 3. Ratios of observed pn+MOS data obtained on 6 March 
2006 (top panel) and 15 October 2006 (bottom panel) and the 
best-fitting absorbed power-law models. 



off is common among bright ULXs and is well described by 



thick plasma. Ja 


Erawal & Misra 20061; Stobbart et al.l 20061; 


iGoad et all 12006 


; iDewangan et al. 2006; Roberts 20071). To 



test such a scenario, we replaced the cutoff PL model with 
the thermal Compton izatio n model (nthcomp ) described by 
IZdziarski etTall (| 19961) and lZvcki et all i|l999l) . The free pa- 
rameters of the nthcomp model are the asymptotic power- 
law index (T t hcom P ), electron temperature (kT e ) and the 
seed photon temperature (kTs). The electron scattering op- 
tical depth (r) can be calculated from the asymptotic power- 
law photon index (T C ompth) and electron temperature (kT e ) 
as follows 



9 



1 



1 1/2 



4 1 (kT e /m e c 2 )r(l + r/3)J " 111 

with Ythcomp = ct + 1 l|Sunvaev fc Titarchuklll980l ). Since a 
soft X-ray excess emission, described by an MCD, is clearly 
detected, it is simplistic to assume that this component 
provides the seed photons for the Comptonization. There- 
fore, the temperature of the soft excess and the seed pho- 
tons were kept the same and varied together. The shape 
of the soft component was assumed to be that of a MCD. 
The MCD+nthcomp model provided a good fit (x 2 /dof = 
717.9/671) and the best-fit parameters are reported in Table 
To confirm that the results do not depend on the Comp- 
tonization model used, we have fitted the data using the 
XSPEC model compTT (iTitarchukll 19941 ; iHua fc Titarchukl 
ll995l ; lTitarchuk fc Lvubarskiill 19951 ) instead of nthcomp and 
obtained a \ 2 /dof = 716.3/671 with similar parameter val- 
ues. In Fig. 21 we show the unfolded EPIC-pn spectrum, 
the best-fitting absorbed MCD+ncompth model in panel (a) 
(upper curve), and the deviations of the data from the model 
in panel (c) . For clarity we show only the EPIC-pn data. Un- 
like XPSEC, the unfolded spectrum in ISIS is derived in a 
model-independent way as follows: 




fvnfold(I) — 



[C(I) - B(I)]/At 
j R(I,E)A(E)dE' 



(2) 



1 2 

Energy (keV) 

Figure 4. Results of spectral fitting to the EPIC-pn data ob- 
tained on high (6 March 2006) and low (15 October 2006) states, 
(a) Lower plot: Observed data for the high state and the best-fit 
model - thermal Comptonization in a cool, optically thick medium 
Unfolded (nthcomp) ; Upper plot: Observed data for the low state 
and the best-fit model MCD+nthcomp. (b) Deviations of the ob- 
served data of the high state from the best-fit absorbed nthcomp 
model, (c) Deviations of the observed data for the low state from 
the best-fit absorbed MCD+nthcomp model, (d) Difference in the 
shape of the spectra for the high and low states. The plot shows 
the deviations of the data of the low state from the best-fit model 
to the high state data but renormalized to match the fluxes. 



where C(I) is the number of total counts in the energy bin 
I, B(I) is the number of background counts, At is exposure 
time, R(I, E) is the normalized response matrix and A(E) 
is the effective area at energy E. This definition produces a 
spectrum that is independent of the fitted model. The best- 
fit model has been over-plotted at the internal resolution 
of the ancillary response function. The spectrum and the 
model are expected to match only i n the particular case of 
unit response matrix l|Nowakll2005l ). Hence there is devia- 
tion between the model and the data points at low ener- 
gies. Note that the model independent unfolded data points 
clearly shows a curvature at high energies, which signifies 
that the curvature is real and not a modeling artifact. 

We then proceed to the spectral analysis of the pn and 
MOS data of the high state. As mentioned earlier in sec- 
tion [2] the extension of the ULX in the EPIC-pn image 
overlaps with the chip gap, hence a good fraction of the 
events at high energies were discarded as these events have 
poorly determined patterns. Consequently, the EPIC-pn and 
the MOS data have comparable signal-to-noise. Simple ab- 
sorbed power-law model fitted separately to the pn, MOS1 
and MOS2 data resulted in similar photon indices and we 
did not find significant differences in the residuals, there- 
fore we present results based on joint spectral fitting of the 
pn and MOS data as before. The simple power law model 
resulted in a minimum \ 2 ~ 455.3 for 440 dof. With a 
reduced x 2 — 1.03, a simple power-law model seems ade- 
quate and does not seem to warrant exploration of more 



6 G. C. Dewangan et al. 



Table 2. Best-fit spectral model parameters for NGC 1313 X-l. 



Parameters 


6 March 2006 


15 October 2006 




nThComp 


MCD+nThComp 




^ . Jul iixeu i 


2 ™+°-i 6 


kT i71 ( eV) a 




227 


JMCD 




1 94+0-15 v in— 12 
0.12 x lu 




122 ± 15 


227_ 12 


kT e { kcV) d 


1 fi7+ - 28 


9 97 +0.18 
z - z '-0.15 


r e 


9 n9 +0.05 
z - u -0.04 


1.67 ±0.03 


7 comp 


7.5lig;« x io- 12 


(3.1 ±0.1) x 10~i 2 


19 
■IX 


7.5ltg;S x IO"" 


4.34^;^ x 10-12 


rh 


17.9 x 10 39 


11.9 x 10 39 


X 2 /ciof 


421.4/439 


717.9/671 



"Inner disk temperature of disk blackbody model which is tied 
to kT s . 

^0.3 — 10 keV flux of the MCD component in ergs cm" 2 s~ 

c Temperature of the seed photon source of the Comptonization 

model 

d Hot electron plasma temperature. 

e Photon Index of the Comptonization model. 

^0.3—10 keV flux of the Comptonization model in ergs cm - 2 s - 1 . 

9 Unabsorbed 0.3 — 10 keV flux in ergs cm - 2 s - 1 . 

^Unabsorbed bolometric luminosity in ergs s - 1 . 



complex models. Indeed iFeng fc KaaretJ (|2006T > were able to 
model several similar quality data of this source with a sim- 
ple power-law. However, the analysis of better quality data 
described above suggests that the spectrum of the source 
could be more complex. Moreover, as we show below, there 
is evidence for a more complex spectral shape even for this 
statistically poorer data set. We have plotted the ratio of the 
pn+MOSl±MOS2 data and the best-fit PL model in Fig- 
ure^ upper panel). The absorbed power-law model is a good 
description of the data below ~ 6 keV. There is a likely cut- 
off or curvature at high energies. To verify the possible high 
energy turnover of the X-ray spectra, we replaced the PL 
component by a cut-off PL. The absorbed cut-off PL model 
yielded A X 2 = -19.9 for one additional parameter as com- 
pared to the absorbed PL model. This is an improvement at 
a significance level of > 99.99% based on the Maximum Like- 
lihood Ratio test. The absorbed nthcomp model provided a 
similar fit with x 2 /dof = 421.4/439. In Figure @] we show 
the unfolded EPIC-pn spectrum, the best-fitting absorbed 
ncompth model in panel (a) marked as 6 March 2006 and the 
deviations in panel (b). Note that, similar to the low state 
data set, the model independent unfolded data points clearly 
show a curvature at high energies which again signifies that 
the curvature is real and not a modeling artifact. It is also 
clear from the model independent unfolded data that both 
the flux and spectral shape varied between the two data sets. 
This is highlighted in panel (d) of the Figure which shows 
the deviations of the low state data from the renormalized 
best fit model of the high state. We have also investigated 



the possible presence of a disk component similar to that 
observed in the low state. Addition of an MCD component 
to the nthcomp model with MCD kTi„ tied to the nthcomp 
kTs did not improve the fit. We calculated an upper-limit of 
Lmcd < 1.4 x 10 39 ergs s -1 for the bolometric luminosity of 
the MCD component in the 0.001 - 100 keV band. We also 
fixed the MCD kTi„ at the value obtained for the low state 
and calculated the upper limit, Lmcd < 6.3 x 10 39 ergs s . 



3 ACCRETION DISK GEOMETRY 

As revealed in Tableful the low state spectrum of the source 
can be described with a thermal Comptonization and a disk 
black body model. The luminosity of the disk black body 
(i.e. the MCD model) is Lmcd ~ 4x 10 39 ergs/s, while that 
of the thermal Comptonization is ~ 7.7 x 10 39 ergs/s. Using 
the XSPEC model function for the model "nthcomp" , one 
can estimate the amplification factor, A = L c / L s ~ 4.8, 
where L c and L a are the luminosities of the Comptoniza- 
tion component and the seed photon input. This allows an 
estimate of L s ~ 1.6 x 10 39 ergs/s. This is consistent with 
an accretion disk geometry where a standard cold accre- 
tion disk is terminated at radius Rt r and the inner region 
is an hot plasma. The inner hot plasma Comptonizes pho- 
tons from the outer disk producing high energy photons. For 
consistency, the seed photon spectral shape for the Comp- 
tonization model is taken to be a disk black body emission 
at the same temperature as the disk emission observed. In 
this model, the fraction of the cold disk photons that enter 
the inner hot region is related to the solid angle which the 
hot region makes with the outer disk. The observed fraction 
of the luminosity of the seed photon component as compared 
to the disk luminosity, Ls / Lmcd ~ 0.4 suggests a solid an- 
gle of Afi ~ 0.4 x 2tt ~ 0.8tt. Note that Ls/Lmcd £ 1 is 
a consistency check for the geometry proposed. In this ge- 
ometry, it would be difficult to reconcile Ls > Lmcd and 
Ls « Lmcd would require an unphysical small subtend- 
ing solid angle. 

The normalization of the disk component provides an 
estimate of the transition radius 
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where i is the inclination angle and k is the color factor. 
In th e standard accretion disk theory (Sha kura fc Sunvaevl 
I1973T ). the luminosity of a truncated disk is given by 
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where M is the mass of the black hole and M is the ac- 
cretion rate. The total bolometric luminosity, Lt = {Lc — 
Ls) + Lmcd can be expressed as Lt = nMc 2 where r\ is the 
efficiency of the system. Thus, the mass of the black hole 
can be estimated to be 



2c 2 Lmcd 
~ 200 M (JL) 
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The Eddington luminosity for a 200M Q solar mass black 
hole is LEdd = 2.5 x 10 40 ergs/s and hence the source is 
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radiating at an Eddington ratio of Lr/L-Edd ~ 0.37. The 
Schwarzschild radius turns out to be r s = 2GM/c 2 ~ 6 x 10 7 
cm, making the transition radius Rt r ~ 17r s . 

Thus the spectral shape of the low state is totally con- 
sistent with a scenario where a ~ 2OOM0 solar mass black 
hole is surrounded by a standard accretion disk which is 
truncated at ~ 17r s . The disk produces a multi tempera- 
ture emission which is observed as a soft component. The 
inner hot region is geometrically thick and subtends a solid 
of angle of ~ 0.8-7T to the disk outside. This allows it to 
Comptonize photons from the outer disk producing the main 
Comptonized component. The system is radiating at a Ed- 
dington ratio ~ 0.4 with a standard radiation efficiency 
77 ~ O.f. 

The spectrum of the high state can be adequately fit- 
ted with a single thermal Comptonization model. The lu- 
minosity of the model is ~ 16.6 x fO 39 ergs/s and an Am- 
plification factor, A ~ 2.5. The input seed photon lumi- 
nosity turns out to be ~ 6.5 x fO 40 ergs/s. This is signif- 
icantly higher by a factor of ~ 4 of the seed photon lu- 
minosity for the low state. If, like in the low state, the 
seed photons are due to an external truncated accretion 
disk, the direct emission from such a disk is expected to 
be Lmcd ~ ~ f .6 x I0 41 ergs s~ 4 which would have 
been easily detected. The upper-limit to the observed disk 
emission, however, is only 6.3 x 10 39 ergs s . Thus the ab- 
sence of a soft component in the spectrum argues against 
such a geometry. Moreover, if the input seed photon shape 
is taken as disk black body emission, the inner disk temper- 
ature turns out to be 0.2 keV which is similar to that of the 
low state, despite the factor of four change in luminosity. 
The luminosity of a truncated accretion disk is oc M / R tr 
while the temperature is oc Af 1/4 /i? 3 { 4 . Thus a near con- 
stant temperature accompanied by a factor of four increase 
in luminosity would imply that the accretion rate M has 
increased by a factor of ~ 8 while the truncation radius has 
increased by a factor of ~ 2. The large factor of ~ 8 change 
in the accretion rate is not reflected in the change in the to- 
tal luminosity of the source, which is only a factor of ~ 1.5 
higher than the low state. Thus a truncated accretion disk 
with an inner hot region does not seem to be a viable model 
for the high state of the source. 

An alternate geometry for the high state spectrum is 
that of a corona covering an optically thick and geomet- 
ricall y thin disk l|Haardt fc Maraschil Il993l : iLiang fc Pried 
1977). In this sandwich model, a fraction of the total gravi- 
tational power is dissipated in the disk, Pd and the rest in the 
corona P c . The disk photons get Comptonized by the overly- 
ing corona and hence only a thermal Comptonized spectrum 
is observed. A fraction £ < 0.5 of the Comptonized photons 
impinge back on to the underlying disk and are absorbed. 
The luminosity of the seed photon is then L s = Pd + £L C , 
where L c is the luminosity of the Comptonized photon. The 
observed total luminosity L — L c (l— £) = Pd + P c - Defining 
the Compton amplification factor A = L c /L s , the dissipated 
powers can be related to the total luminosity as, 

^ = jr~^A Lo < 6 > 

and 

P. = ^jL (7) 



The geometry is valid only for £ < (, ma x = 1/A. In the 
absence of feedback, f = 0, Pd = L /A = Ld and P c — 
L ~Pd- In the other extreme of maximal feedback £ = t; ma x , 
P d = and P c = L . 

For the high state, the spectrum can be interpreted in 
terms of the corona model, with the seed photon luminosity 
assumed to be a black body. The seed photon temperature 
turns out to ~ 0.12 keV and an amplification factor of A ~ 
2.54 which imposes an upper limit on £ mM ~ 0.4. Thus with 
a £ < 0.5 and absence of any other strong soft component, 
the spectrum is consistent with the corona model. The size 
of the corona region R c can be estimated using 
L s = L /(l - i)A. For L ~ 1.7e40 ergs/sec, £ ~ 0.25 and 
A ~ 2.5, R c ~ 3x 10 9 cm. For a 2OOM solar mass black hole 
this corresponds to ~ 50r s . If there is an standard accretion 
disk beyond this radius, its temperature would be < 0.1 
keV and would have a luminosity < 2 x 10 39 ergs/s. Such 
a emission would not have been detected. The Eddington 
ratio for this state turns out be ~ 0.7 as compared to ~ 0.4 
for the low one. 

The temporal variability of the high state spectrum al- 
lows us to further constrain and verify the corona geome- 
try. The fractional root mean squared (rms) variability as 
a function of energy is shown in Figure {5). The simplest 
interpretation would be that the observed variability is in- 
duced by variations in the coronal power P c . From equations 
(0 and 0, the amplification can be written as 

A(T e ,T 3 ,r)= p d + P ; (8) 

rd + C^c 

and the seed photon luminosity as 

T Pd ^>P° rj^A / n \ 

Ls = 1 - <x T s (9) 

A variability in P c (t) = -P c o(l + SP c (t)), will induce cor- 
responding variation in the amplification, 8 A and the seed 
photon temperature ST S . Assuming that the optical depth r 
of the corona does not vary, the above two equations deter- 
mine the time variability of T e (t) and T„(t) as a function of 
8P c (t). The time dependent spectrum can then be obtained 
from T e (t) and T s (t). The standard deviation of the time de- 
pendent spectrum at different energy bands then provides 
a predicted rms variation that can be compared with the 
observed values. For £ = 0.25, (solid line in Figure 0, the 
predicted curve is consistent with the observations. For com- 
parison, curves for £ = (dashed line) and £ = £, m ax = 0.39 
(long dashed line) are also plotted in the figure. Note that 
apart from the overall normalization, the predicated rms 
versus energy plot depends only on £, since the rest of the 
parameters are determined by the time-averaged spectral 
fitting. Thus the coronal geometry naturally explains the 
energy dependent variability of the source. 

As the coronal power increases, the electron tempera- 
ture should also increase making the spectrum harder. Thus 
the model predicts a correlation between flux and hardness 
ratio which is indeed seen in the variability. In the left panel 
of Figure ([1} the hardness ratio is observed to increase with 
flux. While the qualitative behavior is consistent, a quantita- 
tive evaluation is unfortunately not possible. Such a quanti- 
tative estimate would require a fairly accurate knowledge of 
the response of the instrument to convert the model photons 
into counts. As mentioned earlier, for this observation, the 
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Figure 5. The observed rms variation with energy for the low 
state data. The predicted rms variation for £ = 0.25 (solid line), 
£ = (dashed line) and for £ = t; max = 0.33 (long dashed line) 

source extension overlaps with the EPIC-pn chip gap and a 
large fraction of events have uncertain pattern. These events 
with poor spectral calibration although excluded from the 
spectral analysis were included in the temporal analysis. Ex- 
cluding them from the temporal analysis renders the statis- 
tics too poor but their inclusion means that their energies 
are not well constrained. 

The spectrum and rms variability of the source in the 
high state is totally consistent with a geometry where a 
corona covers an underlying disk in a sandwich model. The 
photons emitted from the disk are Comptonized by the 
corona and only a single thermal component is seen. A frac- 
tion £ ~ 0.25 of the Comptonized photons impinge back on 
the disk and are absorbed. The size of the coronal region is 
~ 50r s for a 200Mq black hole and the Eddington ratio is 
~ 0.7. 



4 SUMMARY AND DISCUSSION 

We have investigated the spectral and temporal behavior of 
the bright ULX NGC 1313 X-l using two long XMM-Newton 
observations. On 6 March 2006, the ULX was in a high flux 
state with an average count rate of ~ 1.5 counts/s, while 
on in October 2006 the count rate had decreased to ~ 0.78 
counts/s. The main results of the analysis are as follows. 

The spectrum of the low state shows definite evidence for 
a high energy curvature. It can be fitted with a thermal 
Comptonization model with an optical depth r ~ 15 and 
temperature kT e ~ 2.3 keV. An additional soft component, 
which can be modeled as a disk black body emission with a 
inner disk temperature of kTi„ ~ 0.23 keV, is required for 
the spectral fit. The flux from the source is steady with an 
upper limit on the rms variability of < 3%. 

The spectral shape during the high state is distinctively 



different from that of the low state. The entire spectrum can 
be fitted by a single thermal Comptonization model with an 
optical depth r ~ 13.5 and temperature kT e ~ 1.67 keV. 
There is no evidence for any additional soft component. The 
source is highly variable in this state with an rms variabil- 
ity of 14%. The associated temporal and spectral difference 
between the two observations indicates that the source un- 
derwent a true state transition. 

For a standard radiative efficiency of ~ 0.1, the spectrum 
of the low state is completely consistent with a model where 
a ~ 2OOM0 black hole is surrounded by a standard accretion 
disk truncated at a radius of ~ 17r s (r s = 2GM/c 2 ). This 
outer disk gives rise to the soft component. There is an inner 
hot region which subtends a solid angle of ~ 0.87T to the 
outer disk. Hence photons from the outer disk enter the 
inner hot region, get Comptonized, and are observed as the 
primary thermal Comptonized component. The system has 
an Eddington ratio of Lt /LEdd ~ 0.4. 

The spectrum of the high state is not compatible with a 
truncated disk and inner hot region geometry. This is pri- 
marily inferred by the absence of any strong soft component, 
despite a four fold increase in the input soft photon flux at 
nearly the same temperature. The difference in disk geome- 
try further confirms that a true state transition has occurred 
between the two observations. 

The spectrum as well as the rms variability as a function 
of energy of the high state, is completely consistent with 
a model where a geometrically thin disk has a hot corona 
overlying it. Photons from the underlying disk get Comp- 
tonized in the corona and a single thermal Comptonization 
model is observed. A fraction ~ 0.25 of the Comptonized 
photons impinge back on the disk and are absorbed. The 
observed variability can be naturally ascribed to variations 
in the coronal power. The size of the coronal region is ~ 50r s 
and the Eddington ratio during this state is ~ 0.7. 

The low state of NGC 1313 X-l seems to be analogous 
to the low hard state of Galactic X-ray binaries, which can 
also be broadly explained in terms of a truncated accretion 
disk and an hot inner region. Similarly the coronal geometry 
for the high state is the same as the one invoked for the high 
soft state of Galactic X-ray binaries. However, this is an im- 
portant difference between the hot region/corona of Galactic 
sources and NGC 1313 X-l. The hot region/corona of Galac- 
tic sources is optically thin (r ~ 1) and has a high tempera- 
ture, ~ 50 keV. The hot region/corona for NGC 1313 X-l, is 
optically thick (r > 10) and has a much lower temperature 
~ 2 keV. This difference is more than just quantitative. Op- 
tically thin plasmas with electron temperatures ~ 50 keV 
can be in a two-temperature state, with the proton temper- 
ature being muc h higher (a factor o f ~ 1000) than that of 
the electron gas l|Shapiro et al.|[T976l ). Such plasmas are gas 
pressure dominated. Optically thick and low temperature 
plasmas cannot have such a two-temperature structure and 
would be radiation pressure dominated. Thus, such plasmas 
would have to be dynamically different from those observed 
in Galactic sources. 

With an Eddington ratio of ~ 0.5, the only difference 
between NGC 1313 X-l and Galactic black hole systems 
seems to be the larger mass of the black hole. However, 
the X-ray spectra of bright AGN reveal the presence of an 
optically thin and hot plasma similar to those inferred for 
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Galactic sources. It is not clear why AGN with a ~ 1O 6_9 M0 
black hole mass have similar hot region/corona as ~ IOMq 
Galactic black holes, but ULX with intermediate masses 
have qualitatively different plasmas. The broadband X-ray 
spectra of BHBs in their strongly Comptonized, very high 
state (VHS) require thermal Comptonization in an opti- 
cally thick (r ~ 2 — 3) and relatively cool (kT e ~ 10 keV) 
plasma in addition to multi-color disk blackbody emission 
from a truncated disk and a p ower-law extending to high en- 
ergies l|Done fc Kubotall200rJ ). The curvature produced by 
the low temperature thermal Comptonization component in 
BHBs could be similar to the ~ 5 keV curvature seen in 
the ULX spectra. Furthermore, many narrow-line Seyfert 1 
galaxies and quasars with high Eddington ratio show soft X- 
ray excess emission with kT ~ 200 eV above a hard power- 
law. The soft excess emission can be modeled with Thermal 
Comp tonization in an optically thick and cool pl asma (see 
e.g-. lDewangan et al.ll2007l ; iMiddleton et alj|2008h . Perhaps 
the ~ 5 keV curvature observed in ULX, the ~ 0.2 keV soft 
X-ray excess seen in many AGN, and the ~ 10 keV thermal 
Comptonization seen in strongly Comptonized VHS of BHB 
are all physically the same continuum components. There 
may be an additional hard X-ray (> 10 keV) power-law in 
ULX which would then be analogous to the X-ray power- 
law seen in AGN at £ 2 keV and in VHS of BHBs seen at 
^ 10 keV. However, at this stage such a conjecture would 
be speculative. 

An alternate to the interpretation presented here is that 
the source har bors a ~ 5Mq black hole w i th super Ed ding- 
ton accretion (jShakura fc Sunvaevl Il973l ; iKinef 120081 ) . For 
such high accretion rates, the luminosity does not scale with 
M, and instead is given by L ~ Lscid(l + log(M /Me))- 
Moreover, the characteristic radius at which the luminos- 
ity is released is ~ (27 /4)(M / Me)t s , which is much larger 
than the standard value of ~ 5r s . The luminosity is also ex- 
pected to be further enhanced by some moderate beaming. 
Such an interpretation is attractive because (a) it does not 
require exotic formation and accretion mechanism needed 
to explain intermediate mass black holes, (b) the high ac- 
cretion rates required are perhaps expected when there is 
thermal time scale mass transfer from donor stars and (c) 
it explains naturally why the spectra of ULX are distinct 
from Galactic sources and AGN which are sub-Eddington 
sources. However, the absence of a simple representation of 
the expected mo del spectrum does not allow for detailed 
spectral analysis. iKind l|2009r ) associate the soft excess seen 
in some ULX as emission from such a super Eddington disk. 
This is reasonable if the soft excess is the energetically dom- 
inant component and not for cases like the low state of NGC 
1313 X-l, where the soft component is less luminous than 
the main Comptonization one. Another issue is the factor of 
~ 2 luminosity variability between the low and high state of 
NGC 1313 X-l. Since luminosity varies as log(M /Me) this 
would would imply a large accretion rate variation which 
in turn would imply a large variation in the characteristic 
radius oc M. Such a large variation in radius would result in 
a significant change in the characteristi c temperature of the 
disk, which is not observed. King (2009) invoke variations in 
the moderate beaming factor to explain such observed vari- 
ations, which although plausible is rather ad hoc. Finally, 
if the super Eddington accretion model is correct, then the 
self consistent standard disk modeling which explains both 



the spectra as well as the variability versus energy, has to 
be considered as a coincidence. What is required is a better 
theoretical prediction of the super Eddington disk spectra 
which may have to include the possibility of a additional 
dominant Comptonization component. 

In the standard disk framework, there is also a need 
for better theoretical understanding of a radiation pressure 
dominated corona as well as a hot inner region. The dynam- 
ics of such structures needs to be studied and ascertained 
that they are viable stable systems. Perhaps more impor- 
tantly, the results presented in this work need to be con- 
firmed by longer observations of the same source as well as 
other sources which show similar characteristics. The latter 
is important since there may be different kinds of ULX, even 
among the brightest ones. 
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